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The electric quadrupole moment for the 5d^D3/2 state of ^^^Yb+j has been calculated using the 
relativistic coupled-cluster method. Earlier a similar calculation was performed for the 4d -D5/2 
state of ^^Sr"*" which is the most accurate determination to date [PRL, 96, 193001 (2006)]. The 
present calculation of the electric quadrupole moment of ^'^^Yb"'" yielded a value 2.157eao where the 
experimental value is 2.08(ll)eao; 0,0 is the Bohr radius and e the elementary charge. We discuss in 
this paper our results in detail for ^^^Yb+ and highlight the dominant correlation effects present. We 
have presented the effect of inner core excitations and their contribution to the electric quadrupole 
moment, which is a property sensitive to regions away from the nucleus. 

PACS numbers: 31.15.Ar, 31.15.Dv, 32.30.Jc 



I. INTRODUCTION 

The frequency of any periodic event like the mechan- 
ical oscillation of pendulum or stable atomic frequencies 
can be used to define the unit of time. The frequencies 
derived from selected atomic resonant transitions are par- 
ticularly preferred due to various advantages they offer 
compared to mechanical oscillations. They are extremely 
stable, accurately measurable and reproducible. Though 
the cesium atomic clock frequency [l|, 01 is accurate to 4 
parts in 10^^, a variety of atoms and ions have been pro- 
posed as candidates for the next generation of frequency 
standards 0, Of these candidates, certain ions with 
forbidden transitions in the optical regime are of special 
importance. Trapped and laser cooled single ions and 
neutral atoms trapped in optical lattices are currently 
the leading candidates for atomic clocks [1, 0| . Frequency 
standard experiments with trapped ions require consid- 
erable skill and ingenuity. It is a indeed very challenging 
to measure the clock frequencies to a high degree of pre- 
cision. Effects like the second-order Zeeman, the electric 
quadrupole shift, etc. arise from the interaction of the 
ion with stray fields. It was shown recently that the sys- 
tematic effects caused by these shifts will not limit the 
accuracy of the optical clock 

It was shown earlier that Yb+ is one of the suitable 
candidates for defining a frequency standard |8|j9| . Other 
candidates like Sr+ fm,^, Ca+ [H, Ba+ [li|7Hg+ Q 
etc have also been considered for setting up the frequency 
standard. In particular, Yb"*" is a very versatile candi- 
date, having the clock transitions in the visible, IR and 
microwave regions. The transition which is being con- 
sidered for the frequency standard experiment is the for- 
bidden electric quadrupole {E2) transition between the 
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ground state {6s^Si/2, F ~ 0) and the metastable ex- 
cited state {5d'^D3/2: F = 2). Precise measurements of 
the electric quadrupole moments of these ions have been 
performed [3, d, [^, [13, Comparison of the experi- 
mental values of electric quadrupole moments with those 
calculated theoretically, would serve as excellent tests 
of relativistic atomic theories. An earlier calculation of 
this quantity for Sr+, using relativistic coupled-cluster 
theory yielded a value (2.94 ± 0.07) eag for the 4d 
^^5/2 state, which was in very good agreement with it's 
measured value (2.6 ± 0.3) [1J|, where e is the electronic 
charge and oq is the Bohr radius. This was the most ac- 
curate calculation of the quantity for the 4d state 
of Sr+ to date and the excellent agreement of the mea- 
sured and the calculated values indicates the potential of 
the method used. In this paper, the electric quadrupole 
moment of Yb"*" in the state bd?D^i2 has been calculated 
using relativistic coupled-cluster (RCC) theory. Indeed 
a comparison of this property of a heavy ion like Yb"*" 
with accurate experimental data is a far more stringent 
test of RCC than the corresponding comparison for Sr+. 
The calculation of EQM for Yb+ is computationally more 
demanding due to the presence of a large number of oc- 
cupied orbitals. In our calculation, the entire core has 
been excited. This leads to a rapid proliferation in the 
number of cluster amplitude equations with the size of 
the virtual space considered, and therefore a very large 
increase in the number of computations necessary to de- 
termine these amplitudes. Obtaining convergence for the 
large number of cluster amplitude equations with an ap- 
propriate iterative method is a daunting task for heavy 
atomic systems like Yb+. The calculations on such sys- 
tems hence involve the combination of the power of the 
relativistic many-body theories with the state-of-the-art 
high performance computational techniques. 

An outline of the application of the RCC method to 
calculate atomic electric quadrupole moments has al- 
ready been presented in The details of this theory 
have been discussed in several papers Here we 
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shall give the salient features of the method for complete- 
ness. This paper is organized as follows : Section ITll and 
Section lUll deal with the theoretical methods we have 
employed and the details of our calculation respectively. 
In Section IIVI we present our results and compare with 
the available data. We have also discussed the effects of 
different many-body contributions. Finally we conclude 
in Section [V] and highlight the important findings of our 
work. 



II. THEORETICAL METHODS 

A. Relativistic coupled-cluster theory for closed 
shell atoms 

We start with the iV-electron closed-shell Dirac-Fock 
(DF) reference state |$), which is the Fermi vaccuum for 
the present formulation. In coupled-cluster (CC) theory 
the exact wavefunction for the core sector in terms of this 
reference state is given by, 

|*> -exp(r)|$), (1) 

where T is the cluster operator which takes into account 
the excitations from the closed-shell core to the virtual 
orbitals. In singles and doubles (SD) approximation, the 
cluster amplitude T is written as 

T = Ti+T^ - E ^a + \ E H<^^<'<^) €1 (2) 

ap abpq 

Ti and T2 being the cluster amplitudes for single and 
double excitations respectively and the curly brackets 
denote the normal ordering with respect to the Fermi 
vaccuum. This is known as coupled-cluster with sin- 
gles and doubles, namely CCSD. Here and t^l are the 
corresponding single particle amplitudes and a,b..{p, q..) 
stand for occupied (virtual) orbitals and {• ■ • } denotes 
normal ordering with respect to the common reference 
state (vacuum) |$). For a one- valence one-dimensional 
model space the label 'v' is used to represent a valence 
orbital. In our approach we deal with the normal ordered 
Hamiltonian which is defined as 



(EA) we obtain the (iV -t- l)-electron open shell system 
as 



Atom(0, 0) + e — > lon(0, 1) . (4) 

For a single valence system, we start with the reference 
state 

K^')^al\^) (5) 

where v denotes the valence orbital as mentioned in the 
previous section and the operator represents creation 
of a particle in the valence space. The many-body exact 
open-shell wavefunction for the {N + l)-electron open 
shell system now becomes, 

1*^+1) =exp(T){exp(5.)}|<i>f+i), (6) 

where the curly brackets denote the normal ordering with 
respect to |$). For a single valence system, the operator 
exp(S'i,) turns out to be (1 -I- S^) 

|vI/f+i)=exp(T){(l-H5.)}|<i>f+i), (7) 

with 

Sv = Slv + S2v = E Wp°'^} + ^ E WpO'l'^bav} s^l . 
v^p bpq 

(8) 

Here corresponds to the excitation operator in the va- 
lence (v) sector and and and s^' are the singles and 
doubles amplitudes respectively. The evaluation of the 
cluster amplitudes are discussed elsewhere [? ]. Apart 
from singles and doubles, only approximate triple excita- 
tions (CCSD(T)) have been included. In this calculation, 
we have used OSCC-EA to obtain the 5d ^1)3/2 state of 
^^^Yb"*" which is followed by property calculations as given 
in sub-section III CI 



C. Calculation of Expectation values 



■H = H-{^H\^)=H-Edf, (3) 
where Edf is the Dirac-Fock energy. 



B. Open shell coupled-cluster theory for single 
valence system : Electron Attachment (OSCC-EA) 

To determine the wavefunctions for the open shell or- 
bitals we employ open-shell coupled-cluster method for 
electron attachment (OSCC-EA) for the valence particle 
(O/i, \p) sector. Using the scheme of electron attachment 



The expectation value of any operator O with respect 
to the exact state is given by 

(o) = — 

_ ($^+i|{l + gt}o{i + g}|^jv+i) 

~ ($^+1 1 {1 + 5t} exp(rt) exp(T) {1 + S} \<pN+iy 

where O = exp{T^)0 exp{Tc) is the dressed operator. 
The first few terms of the operators in the _above_ ex- 
pression (Eq. Q) can be identified as O, OSi, OS2, 
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SlOSi etc. The corresponding matrix elements are re- 
ferred to as dressed Dirac-Fock (DDF), dressed pair cor- 
relation (DPC) and dressed core polarization (DCP) re- 
spectively. We use the term 'dressed' because the opera- 
tor O includes the effects of certain core excitations, i.e., 
core-correlation effects. In addition to the above, we can 
identify a few other terms which play a non-negligible 
role in determining the correlation effects. One of those 
terms is SlOSi +c.c which we call the dressed higher or- 
der pair correlation (DHOPC) since it directly involves 
the correlation between a pair of electrons. Diagrams 
representing these terms have already been presented in 
ref. but we have nevertheless given them here for 
clarity. 
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Figure 1: The diagrams (a) and (c) are subsets of dressed pair 
correlation (DPC) diagrams. Diagram (b) is one of the direct 
dressed core-polarization (DCP) diagram. 



where the sum is over all the electrons and z is the co- 
ordinate of the jth electron. To calculate the quantity 
we express the quadrupole operator in its single particle 
form as 



(2) 



(11) 



and the single particle reduced matrix element is ex- 
pressed as [1^ 



O/ll lln^ m = O/ll C^i^ m / dr {VfV, + QfQ,) 



(12) 

In Eq. l|12p . the subscripts / and i correspond to the final 
and initial states respectively; P and Q are the radial part 
of the large and small components of the single particle 
Dirac-Fock wavefunctions respectively and ji is the total 
angular momentum for the ith electron. The angular 
factor is given by 



{Jf\\C^U\J^)= (-l)(^/+i/2)y'(2j/ + i)y(2j7TT) 



where 



n{l,k,l') 



-1/2 1/2 )-(^'^'^') (13) 



1 iil + k + l' even 







otherwise 



D. Electric Quadrupole Moment 

The interaction of the atomic quadrupole moment with 
the external electric-field gradient is analogous to the 
interaction of a nuclear quadrupole moment with the 
electric fields generated by the atomic electrons at the 
nucleus. In the presence of the electric field, this gives 
rise to an energy shift by coupling with the gradient of 
the electric field. Thus the treatment of atomic electric 
quadrupole moment is analogous to its nuclear counter- 
part. 

The quadrupole moment of an atomic state 
J, M)) is defined as the diagonal matrix element 
of the quadrupole operator with the maximum value Mj 
and is expressed as 



(10) 



Here 7 specifies the electronic configuration of the atoms 
which distinguishes the initial and final states; J is the 
total angular momentum of the atom and F is the sum- 
mation of nuclear and atomic angular momentum with 
Mp its projection. The electric quadrupole operator in 
terms of the electronic coordinates is given by. 



I and k being the orbital angular momentum and the 
rank respectively. 

Finally using the Wigner Eckart theorem we define the 
electric quadrupole moment in terms of the reduced ma- 
trix elements as 

0/1 ©1^^ m = i-iy^-"-' ( I ) (jfW e(2) m 

(14) 



III. 



COMPUTATIONAL DETAILS 



This calculation is performed in the following steps : 
The first step being the generation of single particle basis 
for Yb"'""'" using the Gaussian basis set expansion. This 
is followed by the generation of the coupled cluster am- 
plitudes (T) for the closed-shell Yb++ system. In the 
next step, the virtual orbitals 6s and 6(^3/2 are generated 
using the open-shell coupled cluster method for electron 
attachment (OSCC-EA). This is followed by the property 
calculations as given in subsection III CI 

The orbitals used in the present work are generated by 
kinetically balanced finite basis set expansion (FSBE) of 
Gaussian type orbitals (GTO) [20] 



F^,k{r) = r'=exp(-air2). 



(15) 
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with fc = 0, 1, 2 • ■ • for s,p, d, • • • type functions, respec- 
tively. The exponents are determined by the even tem- 
pering condition [2l| 

a,^aG(3'-\ (16) 

The starting point of the computation is the generation 
of the Dirac-Fock (DF) orbitals pS] which are defined on 
a radial grid of the form 

r, = ro[exp(i-l)/i-l] (17) 

with the freedom of choosing the parameters tq and h. 
All DF orbitals are generated using a two parameter 
Fermi nuclear distribution 

''^ l + cxp('(r-c)/a)' ^^^^ 

where the parameter c is the half charge radius and a 
is related to skin thickness, defined as the interval of the 
nuclear thickness in which the nuclear charge density falls 
from near one to near zero. Table U contains the infor- 
mation about the basis functions used in the calculation 
to determine the electric quadrupole moment of 5d ^1)3/2 
state of "iYb+. 

Table I: No. of basis (NB) functions used to generate the even 
tempered Dirac-Fock orbitals and the corresponding value of 
ao — a X 10~^ and /3 used. NP and NH stand for number or 
particles and number of holes respectively. 
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IV. RESULTS AND DISCUSSION 

The contribution of the important physical effects to 
the electric quadrupole moment of ^''^Yb'*' is given in 
TableHH The total value that we have obtained for this 
quantity 83/2 (calculated) — 2.157eag is within the er- 
ror bounds of the measured value 83/2 (measured) = 
(2.08 ± 0.11)ea§ 0]. It would be instructive to com- 
pare our present calculation based on the RCC theory 
with a previous calculation performed by Itano using the 
relativistic configuration interaction (RCI) method [2a|. 
The value obtained by Itano is O3/2 = 2.174ea§ |22| . 
The RCI calculation of Itano uses a multiconfiguration 
Dirac-Fock (MCDF) extended optimized level (EOL) or- 
bital basis. The configurations included in the latter 
calculation constitute a subset of the configurations in 
our calculation. In particular, they correspond to the 
correlation effects arising from the single and double ex- 
citations from the core and the valence, i.e., terms in- 
volving Ti, T2, Si and S2 in our RCC calculation. In 



our calculation, the above effects have been included to 
all orders. The virtual orbitals considered by Itano were 
{(6 - 10)s, (6 - 10)p, (6 - 10)d, (5 - 7)/, (5 - 6)g, 6h}. The 
calculation carried out by Itano involves 2 main steps. In 
the first step, the SPOs are obtained by minimizing an 
energy functional in a limited orbital space in the frame- 
work of MCDF-EOL. The second step involves a fairly 
large RCI calculation to account for the correlation ef- 
fects. Itano's calculation incorporates the valence-core 
correlation (single and the double excitations) arising 
from the {5d} and {4f,5s,5p} shells and single excita- 
tions from {4s,4p,4d,3d} (core-core correlations). Inspite 
of these differences, the results are in good agreement. 

For each symmetry we have considered more virtual or- 
bitals than Itano. In addition, we have considered single 
and double excitations from all the core orbitals where 
Itano has considered only single and double excitations 
from 5d and {4/, 5s, bp} core orbitals, but not more than 
one single core excitation at a time. 

From Table. [Ul we see that the dressed Dirac-Fock 
contribution is the largest and it is a substantial fraction 
of the total EQM, inspite of the large number of core- 
valence excitations. The second largest contribution is 
from the DPC effects and third being DCP effects. DPC 
effects arise from the terms like OSi and SlOSi. Sub- 
section lTl Cl gives the details of these different many-body 
terms. 5*1 is an operator of rank and the valence orbital 
in ^''^Yb^ is a 6^3/2 orbital and hence it excites the va- 
lence electron to a virtual orbital of the same parity and 
angular momentum giving a large contribution through 
the virtual d orbitals. Clearly the effects of the valence 
5d3/2 excitations are the most important in the case of 
i7iYi3+_ ig interesting to note that though ^^^Yb"*" 
has more filled shells than for **Sr+, the contribution of 
the DDF, DPC, DCP terms follow the same trend for 
i71yi-,+ and ^^Sr"*". The core-virtual electric quadrupole 
excitations/deexcitations contributing to the DCP dia- 
grams involve excitations from f - f or d - d or f - p, etc 
orbitals. The DCP diagrams involve the contributions 
from the O and the 5*2 matrix elements. Though the O 
matrix element could be large in some cases, the contri- 
bution of the matrix element of the product OS2 turns 
out to be two orders of magnitude smaller than that of 
the DDF term, inspite of the possibility of large number 
of core- virtual excitations. Similar trends were observed 
for the EQM of ^^St+ . Table [l] has more details about 
the virtual orbitals and the active space. From Table, [ill 
we find that the contributions of the terms DDF, DPC, 
DCP and DHOPC to the final value are 119%, -14%, - 
2.2%, 1.2% respectively for ^^St+ and are 116%, -13.3%, 
-1.3%, 1.34 % respectively for ^'^^Yh+ . Also, the total 
contribution of DPC, DCP and DHOPC to DDF is -1.1 
%. DPC is ~ 11% of the DDF. DCP is ~ 1% of DDF for 
the case of ^''^Yb+. It has been observed that this trend 
is similar for ®^Sr+ even though ^^^Yb+ has a larger core. 
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Table II: Contributions from the electric quadrupole moment 
(in eag) of the Sd^Dg/a state of ^''^Yb+ and the 4:d^Dr,/2 
state of **Sr+, corresponding to different many-body effects 
in the CCSD calculation. The terms like DDF, DCP, DPC, 
DHOPC are explained in the text. 



Ion DDF DPC DCP DHOPC Total Expt. [9| 



Yb+ 2.500 -0.287 -0.0280 0.029 2.157 2.08±0.11 
Sr+ 3.496 -0.4306 -0.0642 0.0353 2.94 2.6± 0.3 



nation of EQM of the Sd^Dg/astate of ^'^'^Yh+ is within 
the experimental Hmits. It highhghts the abihty of RCC 
to capture the interplay between the relativistic and the 
correlation effects in heavy single valence ions. We have 
also determined the various leading many-body effects 
arising in this calculation. To our knowledge this calcu- 
lation yields the most accurate theoretical value of EQM 
of the 5(i^il'3/2state of ^^^Yb+to date. It is a useful theo- 
retical supplement to the experimental search for optical 
frequency standards. 



V. CONCLUSIONS 

In summary, we have used relativistic coupled-cluster 
theory to calculate the electric quadrupole moments 
(EQM) of the state of "iYb+. Our determi- 



Acknowledgments 

This work is partially supported by BRNS project 
no. 2002/37/12/BRNS. RKC acknowledges the Depart- 
ment of Science and Technology, India (grant SR/Sl/PC- 
32/2005). 



[1] R. Wynands and S. Wyers, Metrologia, 42, S 64 (2005). 

[2] 'http://tf.nist.gov/cesium/atomichistory.htm' 

[3] H.S. Margolis, et al., Atomic Phys. 20, Proceedings XX, 

ICAP, 2006. Published by AIP Conference Preceedings, 

869. 

[4] L. Hollberg et al, J. Phys. B 38, S469 (2005). 

[5] W. H. Oskay et al, Phys. Rev. Lett., 97, 020801 (2006). 

[6] M. Takamoto, F-L Hong, R. Higashi, H. Katori, Nature 

435, 321 (2005). 
[7] W. H. Oskay, W. M. Itano and J. C. Bergquist, Phys. 

Rev. Lett. 94, 163001 (2005). 
[8] J. Stenger et al, Opt. Lett. 26, 1589 (2001). 
[9] T. Schneider, E. Pein, and C. Tamm, Phys. Rev. Lett. 
94, 230801 (2005). 
[10] J. E. Bernard et al, Phys. Rev. Lett. 82, 3228 (1999). 
[11] H. S. Margolis et al, Phys. Rev. A 67, 032501 (2003). 
[12] C. Champenois et al, Phys. Lett. A 331, 298 (2004). 
[13] J. A. Sherman et al, Phys. Rev. Lett. 94, 243001 (2005). 



[14] G. P. Barwood et al, Phys. Rev. Lett. 93, 133001 (2004). 
[15] P. Dube et. al, Phys. Rev. Lett. 95, 033001 (2005). 
[16] C. Sur et al, Phys. Rev. Lett. 96, 193001 (2006). 
[17] R. F. Bishop, Lecture Notes in Physics, Microscopic 

Quantum Many-Body Theories and their Applications, 

p.l, Eds. J. Navarro and A. Polls, Springer- Verlag-Berlin, 

Heidelberg and New York (1998). 
[18] R. J. Bartlett, Modern Electronic Structure Theory, vol- 

II, p. 1047, Ed. D. R. Yarkony, World Scientific, Singapore 

(1995). 

[19] I. P. Grant, J. Phys. B, 7, 1458 (1974). 

[20] R. K. Chaudhuri, P. K. Panda and B. P. Das, Phys. Rev. 

A, 59, 1187 (1999). 
[21] R. C. Raffenetti and K. Ruedenberg, J. Chem. Phys, 59, 

5978 (1973). 

[22] Wayne M. Itano, Phys. Rev. A, 73, 022510 (2006). 



